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ABSTRACT 

M85 optical transient 2006-1 (M85 OT2006-1) is the most luminous member of the small family of 
V838 Mon-like objects, whose nature is still a mystery. This event took place in the Virgo cluster of galax- 
ies and peaked at an absolute magnitude of Mj w —13. Here we present Hubble Space Telescope images of 
M85 0T 2006-1 and its environment, taken before and after the eruption, along with a spectrum of the host 
galaxy at the transient location. We find that the progenitor of M85 0T 2006-1 was not associated with any 
star forming region. The g and z-band absolute magnitudes of the progenitor were fainter than about —4 and 
—6 mag, respectively. Therefore, we can set a lower limit of ~ 50 Myr on the age of the youngest stars at 
the location of the progenitor that corresponds to a mass of < 7 M©. Previously published line indices sug- 
gest that M85 has a mean stellar age of 1.6 ±0.3 Gyr. If this mean age is representative of the progenitor of 
M85 0T 2006-1, then we can further constrain its mass to be less than 2 M . We compare the energetics and 
mass limit derived for the M85 OT 2006-1 progenitor with those expected from a simple model of violent stel- 
lar mergers. Combined with further modeling, these new clues may ultimately reveal the true nature of these 
puzzling events. 

Subject headings: stars: individual (M85 OT 2006- 1 , V838 Mon, M3 1 RV, V4332 Sgr) 



1. INTRODUCTION 

M85 Optical Transient 2006-1 (M85 OT2006-1; 
J122523.82+181056.2) was discovered on 2006 Jan 6 
by the Lick observatory supernova search team (Filippenko 
et al. 2001 8 ) as a faint, V ~ 19.3 mag transient in the galaxy 
M85 (NGC4382), which is at a distance of 17.8 Mpc (Mei 
et al. 2007). Subsequent spectroscopy, as well as visible 
light and infra-red (IR) photometry, presented in Kulkarni 
et al. (2007), showed that M85 0T 2006-1 has a recession 
velocity of 880 ± 130 km s , and is therefore associated with 
M85. Moreover, we showed that the temporal and spectral 
properties of this object are unlike those of supernovae, 
novae, or luminous blue variables. 

M85 0T 2006-1 peaked at absolute /-band magnitude of 
about —13. The light curve settled into a ~ 60 day plateau, 
followed by a decrease in bolometric luminosity during which 
the black-body emission peak shifted toward near-IR wave- 
lengths. The early spectrum of M85 OT 2006-1, obtained six 
weeks after discovery, resembles that of a ~ 4600 K black 
body, with Ha and Hj3 narrow emission lines (full width at 
half maximum of ~ 350km s -1 ), along with several other 
unidentified emission lines. Spitzer IR observations obtained 
about six months after the discovery revealed a ~ 1000K 
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black body spectral energy distribution (Rau et al. 2007). 

The spectral and temporal properties of this object resemble 
those of M31-RV (discovered by Rich et al. 1989; e.g., Mould 
et al. 1990; Bryan & Royer 1992), V838 Mon (discovered by 
Brown 2002; e.g., Kimeswenger et al. 2002; Bond et al. 2003; 
Corradi & Munari 2007), and possibly the less studied object 
V4332 Sgr (Martini et al. 1999). However, the M85 transient 
is the most luminous member of the V838 Mon class. The 
favored model for this emerging class of V838 Mon-like ob- 
jects (also known as luminous red novae 9 ) is that they are the 
result of stellar mergers (e.g., Soker & Tylenda 2006). How- 
ever, other models have been suggested to explain these ob- 
jects (e.g., Retter & Marom 2003; Lawlor 2005). The nature 
of these events, with their energetics lying between the realms 
of supernovae and novae, remains uncertain. 

In this paper, we present Hubble Space Telescope (HST)- 
Advanced Camera for Surveys (ACS)AVide Field Camera 
(WFC) and Near Infrared Camera and Multi-Object Spec- 
trometer (NICMOS) observations, as well as Palomar 5 m 
spectroscopy, of the environment of M85 0T 2006-1. The 
observations are used to characterize the environment of the 
transient and to set a limit on the mass of the progenitor. 

2. OBSERVATIONS 

M85 was observed using HST/ACS on 2003 as part of the 
HST-ACS Virgo Cluster Survey (Cote et al. 2004). Sub- 
sequently, the transient was observed serendipitously with 
ACS/WFC and NICMOS/NIC2 in 2006 (GO-10515) as a 
follow-up study to Peng et al. (2006). The ACS observa- 
tions on 2006 were obtained 18 days after the discovery of the 
transient. The log of observations, the measured magnitude 
of the M85 transient or the ~ 3<7 upper limit at the OT lo- 
cation, as derived from the HST images taken on 2003 and 
2006, are listed in TableQ] The HST images of the galaxy and 
the transient environment are presented in Figures Q] and [2] 

9 this term was introduced by Kulkarni et al. 2007. 
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FIG. 1. — Left: HST/ACS F814W-band image of the galaxy M85, obtained on 2006 Jan 24 (18 days after the discovery). The transient, which is well 
detected, is marked by a circle. Right: HST/ACS F475W-band image of M85 after subtraction of the best fit Sersic model (using GalFit; Peng et al. 2002). 
The subtracted model parameters are: effective radius 389"; Sersic index 3.0; axis ratio (b/a) 0.765; position angle 29.3°; diskiness —0.064. A different set of 
structural parameters is obtained when analyzing the azimuthally averaged profile (Ferrarese et al. 2006). We note, that the rough galaxy subtracted image allows 
us to show the lack of dust and structure in the neighborhood of the transient. Note that the gray scale level stretch in the left panel is about 5.2 times larger than 
in the right panel. The white band in the images is due to the gap between the ACS CCDs. The slit of the Palomar 5 m telescope spectrum (Fig. [5} passes through 
the transient location and the center of the galaxy. 




FIG. 2.— Zoom in on the environment of M85 OT 2006-1 HST/ACS and NICMOS/NIC2 images. The circle, with radius of 1", marks the position of the 
transient. Note that the F475W + F814W is a sum of the F475W and the F814W images. 
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TABLE 1 

Photometry and limiting magnitude 



Date 


Exposure 


Band 


Mag. a 


Limitin 


g mag. b 




s 






apparent 


absolute 


2003 Feb 01 


750 


F475W{g) 




> 26.9 


> -4.5 


2003 Feb 01 


1120 


FS5QLP(z) 




> 25.1 


> -6.2 


2006 Jan 24 


2204 


F475W(g) 


20.57 






2006 Jan 24 


2224 


F8Uw(i) 


18.62 


> 25.3 


> -6.0 


2006 Feb 28 


500 


F160W (H) 


17.82 


> 21.2 


> -10.1 



a Vega based magnitude corrected for infinite aperture (Sirianni et al. 2005). 
Errors in photometry are about 0.02 mag for the ACS observations, and 
0.05 mag for the NICMOS observations. The NICMOS magnitude is cal- 
ibrated using 2MASS stars in the field of view. b Vega based limiting magni- 
tude as estimated by adding artificial point sources to the images in the neigh- 
borhood of the transient and inspection of the images for the added sources. 
The absolute magnitudes are calculated assuming a distance of 17.8 Mpc to 
M85 (Mei et al. 2007) and Galactic extinction of F B -y = 0.031 (Schlegel, 
Finkbeiner, & Davis 1998; Cardelli, Clayton, & Mathis 1989). Note however 
that distance estimates to M85 range between 14 Mpc (Ciardullo et al. 2002) 
to 18.6 Mpc (Blakeslee et al. 2001). 
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FIG. 3. — Two-dimensional spectrum of M85 and the transient environment 
(+30" offset from the galaxy center along the slit), obtained about one year 
after the transient discovery, covering the Ha wavelength region (marked by 
ellipse). No Ha emission is seen in the vicinity of the transient. The spectrum 
is shown before sky subtraction. 



On 2007 January 20, after the M85 OT 2006-1 faded away, 
we obtained a spectrum at the location of M85OT2006- 
1. The spectrum (Figure [3]l consist of 4 x 300 s exposures 
with the double beam spectrograph mounted on the Palomar 
5 m telescope. We used the 600 lines/mm grating blazed at 
9500 A in the red arm. The 2" slit was centered on the nu- 
cleus of M85, at a position angle of 185deg. The position 
angle was chosen such that the location of the transient will 
be included in the slit. 

On 25 June 2005, M85 was observed by the Spitzer space 
telescope with the Multiband Imaging Photometer for Spitzer 
(MIPS). The 70 micron image, with exposure time of 670 s, 
is shown in Fig. |4] 

3. RESULTS 

V838 Mon-like objects are found in both young regions 
(e.g., V838 Mon; Afsar & Bond 2007) and old stellar popu- 
lations (e.g., M31 RV; Bond & Siegel 2006). M85 0T 2006-1 
took place in an early-type galaxy. Therefore, as we explain 
below, it can be used to set an upper limit on the minimal 
progenitor mass that can produce V838 Mon-like eruptions. 

From the 2-dimensional spectrum, shown in Fig. [3] we can 
set an upper limit on the flux of the Ha emission at the lo- 
cation of M85 OT 2006-1 of < 6 x 10" 14 ergs cm" 2 s _1 , at the 
3.5(7 level (equivalent to luminosity of 1.6 x 10 39 ergs s -1 ). 
This corresponds to a star formation rate smaller than about 
10~ 2 M yr~' (Kennicutt 1998) in a radius of — 100 pc 




FIG. 4. — Spizter/MIPS 70 micron image of M85. The plus sign marks the 
visible-light center of the galaxy, while the circle marks the position of the 
transient. 



around the transient location. For comparison, the Ha lumi- 
nosity of the Orion nebula is about ~ 10 41 ergs s (Haffner 
et al. 2003; assuming a distance of 392 pc; Jeffries 2007). 
Therefore, our observations rule out the presence of a promi- 
nent star forming region in this location. Moreover, based on 
the far-IR flux in the region of the transient, obtained from the 
Spitzer/MIPS 70 micron image shown in Fig.|4j we can set an 
upper limit on the star formation rate in this region to be less 
than 10~ 5 M yr" 1 (Kennicutt 1998). The absence of H II 
regions in M85 rule out the possibility that the progenitor had 
a delay (from birth to outburst) of < 10 Myr (corresponds to 
<; 40 M Q ), which is the typical life time of H II regions (e.g., 
Mayya 1995). 

An independent limit on the age and mass of the progenitor 
can be inferred from the absence of stars brighter than z-band 
absolute magnitude M z = —6.2 in the transient environment. 
From the Lejeune & Schaerer (2001) stellar tracks, we find 
that stars older than 50 Myr (and therefore, more massive than 
7 Mq) do not reach z-band absolute magnitudes brighter than 
M z = —6.2 mag (at the red-supergiant stage). We note that 
the z-band stellar track magnitudes were obtained by interpo- 
lation of the / and /-bands. Therefore, we can set a lower limit 
on the age of the most massive stars in the transient environ- 
ment to be J> 50 Myr, which corresponds to mass < 7 M 
(assuming solar metallicity). Otherwise, we were likely to 
detect individual stars in this region. We note that we can 
limit the extinction in the transient location to A, < 0.8 mag, 
based on the Balmer lines ratio, assuming case-B recombina- 
tion (Kulkarni et al. 2007). 

Ferrarese et al. (2006) reported possible faint wisps and 
patches of dust in M85. Moreover, Schweizer & Seitzer 
(1992) reported that M85 is somewhat bluer than typical SO 
galaxies, therefore possibly younger. This claim is supported 
by Terlevich & Forbes (2002) who estimated the age and 
metallicity based on line indices. They have found a mean 
luminosity-weighted age of 1.6±0.3Gyr and metallicity of 
[Fe/H]= 0.44 and [Mg/Fe]= 0.08. We note that the actual 
mean age is probably higher than that indicated by line in- 
dices given that younger populations have higher weight than 
old population. If the mean age is representative of the pro- 
genitor of M85 0T 2006-1, then we can set a lower limit of 
about 1 Gyr on the age of M85 OT 2006-1 progenitor/s. This 
further suggests that the mass of the progenitor/s is probably 
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below 2 M© (the life time of solar metallicity > 2 M stars 
is < 1 Gyr; Lejeune & Schaerer 2001). This limit is based on 
the mean stellar age of this galaxy. However, stars younger 
than 1 Gyr may be present in this galaxy in relatively small 
numbers. 

4. DISCUSSION 

Although several models exist for V838 Mon-like objects 
(e.g., Soker & Tylenda 2003; Lawlor 2005), in the absence of 
detailed simulations, the nature of these objects remain elu- 
sive. A clue to their origin can be derived from their en- 
vironment, luminosity function and rate. Given that only a 
small number of these objects are known, and they were found 
serendipitously in various searches, the luminosity function 
and rate are not well constrained. However, the fact that at 
least two events were observed in our Galaxy (i.e., V838 Mon 
and V4332 Sgr) in the last ~ 13 years suggests that they have 
a higher rate than SNe. We can set a lower limit on their rate, 
of 0.019yr~' L^ l w , at the 95% confidence level, where Lmw 
is the Milky Way luminosity. 

Now we discuss the implications of our observations for a 
specific model for V838 Mon-like objects. Soker & Tylenda 
(2006) presented a model for violent stellar mergers in which, 
prior to the merger, the spins and orbital frequencies of the 
binary star are losing synchronization due to the Darwin in- 
stability (e.g., Eggleton & Kiseleva-Eggleton 2001). They 
found that for a given primary mass, the maximal energy 
production obtained for a binary mass ratio of ~ 1 /50, is 
~ 2.5 x 10 _3 GMj /R\, where G is the gravitational constant, 
and M\ and R\ are the mass and radius of the primary star. 
Given the upper limit on the progenitor mass, based on the 
mean stellar age in M85, < 2 M©, and assuming a main- 
sequence mass-radius relation, R M 7 , the maximum avail- 
able energy in their model is short by a factor of three in the 
total energy production, as compared to the radiated energy 
of M85 0T 2006-1 in the first two months, - 8 x 10 46 ergs 
(assuming a distance of 17.8 Mpc to M85; Mei et al. 2007). 
Moreover, it is expected that a large fraction of the energy will 
go into lifting the outer region of the star rather than radiated 



away. Furthermore, if the primary is an evolved star, then its 
radius will be larger than the radius of a main sequence star 
with the same mass, and the extracted energy will be even 
smaller. This suggests that either more detailed modeling of 
violent stellar mergers is required, or that this event is not the 
result of a violent stellar merger. Another possible solution is 
that the mass of the progenitor is somewhat larger. A larger 
progenitor mass will still be consistent with our upper limit 
of 7 M which is based on the absence of stars brighter than 
7^—6 mag. For example, according to Soker & Tylenda 
(2006) model, a 7 M progenitor can yield ~ 4 times more 
energy than a 2 M progenitor and may explain the discrep- 
ancy. We note, however, that other kinds of instabilities can 
lead to stellar mergers (e.g., in triple systems) and that the 
above comparison is valid only for the specific case discussed 
by Soker & Tylenda (2006). 

Existing hydrodynamical simulations of the common en- 
velope phase in stellar mergers (and also star + neutron star 
mergers) predict that the total dissipated energy is of the or- 
der of that observed in V838Mon and M85 0T 2006-1 (e.g., 
Taam & Bodenheimer 1989; Terman et al. 1995; Terman & 
Taam 1996). Moreover, simulations of the common envelope 
phase predicts that most of the envelope will be ejected in 
the equatorial plane (e.g., Taam & Ricker 2006). Indeed, in 
Rau et al. (2007) we reported evidence suggesting that the ex- 
pansion of M85 OT 2006-1 is asymmetric. However, more de- 
tailed hydrodynamical simulations of the vast parameter space 
available for stellar mergers are needed in order to understand 
these processes and to test if V838 Mon-like objects are in- 
deed the results of stellar mergers. 

To summarize, we show that, in contrast to V838 Mon, but 
similarly to M31 RV, M85 OT 2006-1 was probably produced 
by members of an old stellar population (> 1 Gyr), and that 
its progenitor/s mass was probably < 2M . These constraints 
narrow down the allowed venue of stellar models for the na- 
ture of this event. 

This work is supported in part by grants from NSF and 
NASA. 
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